Individual complex Dirac eigenvalue distributions from random matrix theory and 
comparison to quenched lattice QCD with a quark chemical potential 
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We analyze how individual eigenvalues of the QCD Dirac operator at nonzero quark chemical 
potential are distributed in the complex plane. Exact and approximate analytical results for both 
quenched and unquenched distributions are derived from non-Hermitian random matrix theory. 
When comparing these to quenched lattice QCD spectra close to the origin, excellent agreement is 
found for zero and nonzero topology at several values of the quark chemical potential. Our analytical 
results are also applicable to other physical systems in the same symmetry class. 

PACS numbers: 12.38.Gc, 02.10.Yn 



Introduction. Hcrmitian random matrix theory 
(RMT), which describes systems with real spectra, enjoys 
many apphcations in physics and beyond. Dropping the 
Hermiticity constraint results in matrices whose eigenval- 
ues are, in general, complex. Examples are the Ginibre 
ensembles [I] or their chiral counterparts [2]. Although 
these ensembles describe non-Hermitian operators, they 
have found many applications (see [3] for a recent re- 
view), ranging from dissipation in quantum maps jl] 
over quantum chromodynamics (QCD) at nonzero quark 
chemical potential 5 to the brain auditory response de- 
scribed by nonsymmetric correlation matrices 6 . 

Observables that are typically computed in RMT are 
spectral correlation functions. Alternatively, one can 
study the distributions of individual eigenvalues, pro- 
vided that the latter can be ordered. For RMT with real 
eigenvalues, all such distributions are known and have 
found a variety of important applications. For exam- 
ple, the largest eigenvalue follows the Tracy- Widom dis- 
tribution [7] and appears in the longest increasing sub- 
sequence of partitions ^ or growth processes pj. The 
smallest eigenvalue distribution in chiral RMT has be- 
come a standard tool in lattice QCD to extract the low- 
energy constant (LEG) S that appears in chiral pertur- 
bation theory (chPT) and is related to the chiral conden- 
sate [TU' . This distribution is also sensitive to the gauge 
field topology and can be used to distinguish different 
patterns of chiral symmetry breaking |llj . 

In this paper, we generalize some of these results to the 
case of non-Hermitian chiral RMT in the unitary symme- 
try class. We study the distributions of individual eigen- 
values in the complex plane and derive analytical results 
for the chiral RMT introduced in Ref. [12]. Our main 
focus will be on QGD, but our findings are also relevant 
for other systems with complex eigenvalues in the same 
symmetry class. 

In QGD, a nonzero quark chemical potential /i leads 
to a complex spectrum of the Dirac operator. In the 
large- volume limit, chiral RMT is equivalent fT5] to the 



chiral effective theory for the epsilon- regime of QGD [14], 
which is a particular low-energy limit of the full theory. 
Here, a virtue of ^ is that /i couples to the second 
LEG in leading order of chPT, F, which is related to the 
pion decay constant [15 . A comparison of lattice QGD 
data to individual complex Dirac eigenvalue distributions 
from RMT thus allows us to determine both S and F (for 
related methods, see Refs. [TBI fT7] ). 

Unfortunately, lattice QGD with dynamical fermions 
at 7^ faces a serious difficulty due to the loss of re- 
ality of the action. It is very hard to obtain significant 
statistics in unquenched simulations, and therefore we 
will only compare to quenched simulations below. How- 
ever, for /i < m7r/2 or fi^F'^V < 1 (where is the 
pion mass and V is the volume) the sign problem is not 
severe [XBs, and our method can be used to determine 
F from such unquenched lattice data. Therefore we also 
derive RMT results for unquenched QGD, thus adding to 
the predictions for spectral densities [19] and the average 
phase factor [15] . 

What is known from RMT for individual eigenvalue 
correlations in the complex plane? For the non-chiral, 
unitary Ginibre ensemble the repulsion (or spacing dis- 
tribution) of complex levels was computed in ^ and suc- 
cessfully compared to lattice QGD data in the bulk of the 
spectrum [20] . For maximal non-Hermiticity, the distri- 
bution of the largest eigenvalue with respect to radial 
ordering is also known [3T]. However, in QGD it is the 
eigenvalues closest to the origin that carry information 
about topology and LEGs, and therefore we concentrate 
on these in the following. 

The complex spectral correlation functions of the QGD 
Dirac operator at /z ^ were computed from different 
(but equivalent) chiral RMTs in Refs. [H [111 HI and 
compared to quenched lattice QGD in Refs. [TT] |23]. 
Later, a Dirac operator with exact chiral symmetry at 
/i ^ was constructed [24| ES] and tested against chiral 
RMT for topological charge ly — 0, 1. Here, we compare 
the data of Ref. [23^ to our newly derived individual 
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complex eigenvalue distributions, resulting in a much 
improved signal. For a recent review we refer to Ref. |26j . 

Complex eigenvalue distributions. We start by 
defining the gap probability and the distribution of an 
individual eigenvalue in the complex plane. Suppose the 
partition function Z can be written in terms of N com- 
plex eigenvalues Zj of some operator, with a joint proba- 
bility distribution function (jpdf) 7^({z}), symmetric in 
all its arguments, to be specified. (For simplicity, we 
consider only jpdf's with additional symmetry z ^ — z, 
restricting ourselves to the upper half- plane C+.) The 
complex eigenvalue density correlation functions are de- 
fined as 



Results from RMT. The above considerations hold 
for any jpdf, including the jpdf appearing in the lattice 
QCD partition function in terms of complex Dirac op- 
erator eigenvalues and the jpdf of chiral RMT. We now 
consider the latter. The RMT for unquenched QCD with 
^ ^ we use here is given by 



(7) 



The Vandermonde, Ajv({^'}) = Vu>,{^'l " ^])^ 

coming 

from the diagonalization of complex matrices of dimen- 
sion A^x {N ^v) (we take > for convenience), leads to 
a repulsion of eigenvalues. (For the chiral RMTs corre- 
sponding to adjoint or two-color QCD, the Jacobians will 
be different, leading to different patterns of eigenvalue re- 
pulsion, see, e.g., Ref. f2B^.) The weight w depends on 
dynamical quark flavors with masses rrij (/ = !,..., N f) 
and on the number v of exactly zero eigenvalues (corre- 
sponding to the topological charge), 
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The simplest example Ri{z) is just the spectral density. 
The gap probability Ek[J] is defined as the probability 
that there are exactly k eigenvalues inside the set J and 
N — k eigenvalues in its complement J = C+/J, 
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If all Rk are known, the gap probabilities follow as in the 
real case [271. 
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Parameterizing the boundary dJ of J in C+ as z(r) = 
x(t) -I- iy{T), we can define the probability pk{J,T) for 
k — 1 eigenvalues to be inside J, for the eigenvalue Zk — 
z{t) to be on the contour dJ at r, and for N ~ k eigen- 
values to be in the complement J, 
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(Because eigenvalues repel each other in RMT, the prob- 
ability of finding two eigenvalues at z(t) is zero.) An 
ordering on C+ is induced by a family of sets of increas- 
ing area with mutually nonintersecting contours. Via the 
Riemann mapping theorem, this can always be reduced 
to radial ordering. Definitions ^ and Q are related 
through a variational derivative, 

5Ek[J] 
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Employing the expansion (|3]), we can express the pk{J,T) 
through densities. For example, for the first eigenvalue, 

Pi{J, r) = Ri{z{t)) - d^zi R2{zi,z{t)) 

+ ^ CpZ^ d^Z2 R3{Z1,Z2, z{t)) + ... (6) 



where is a modified Bessel function and p, is the chem- 
ical potential in the random matrix model. The first fac- 
tor in Eq. Q originates from the Dirac determinants. 
The non-Gaussian weight function results from an inte- 
gration over angular and auxiliary variables |12j . For 
fi the Zk are back on the imaginary axis. Complex 
RMT yields the following result for the densities [28], 

fe 

Rk{zu---,Zk) = T\w^^f''''>{ze) det Kn{z,,z*), (9) 

^ ^ l<i j<k ■' 

given in terms of the kernel K]\j{zi, z*) of (bi-) orthogonal 
polynomials with respect to the weight of Eq. ^ . In the 
quenched case {Nf = 0), these are given by Laguerre 
polynomials in the complex plane |12j . All unquenched 
density correlations are given explicitly in Ref. |19j . A 
determinental expression follows for the Ek[J] in terms 
of the kernel operator times the characteristic function of 
J. Eq. (|3| is called its Fredholm determinant expansion. 

As mentioned above, in the limit of large volume V, 
RMT is equivalent to QCD in the epsilon-regime |13j . 
In this regime, the chemical potential, the quark masses, 
and the Dirac eigenvalues are rescaled such that the pa- 
rameters a = 2NfL'^ (= VF^fM^), rif = Nmf (= FSm/), 
and ^fc = Nzk (= VSzk) stay finite in the large- 
N (large-!/) limit. In parentheses, we have given the 
scaling of these parameters in terms of the LECs of chPT. 

Quenched case. In the quenched case, the RMT 
result for the microscopic spectral density pi(^) = 
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FIG. 1: Quenched density pi{^) of Eq. ([T0| (left), and 
quenched pi(^) from Eq. ^ (including the first three terms) 
for circular J (right), both for v — and a = 0.174. 



liniAr^oo -Ri(^ = z/N) /N is given by [HJ 



(10) 

where is a modified Bessel function. The rescaled ker- 
nel giving all correlation functions according to Eq. ([9| 
was derived in Refs. [H [TH]. In Fig. [l] we show as an 
example the density pi{^) and the distribution pi(^) of 
the first eigenvalue from Eq. ^ (in which J is chosen 
to be semi-circular and only the first three terms are in- 
cluded). As in the case of real eigenvalues p7], we see 
that the expansion converges rapidly. Higher-order terms 
merely assure that Pi(^) remains zero for large |^|. 



For increasing a, the quenched density Eq. ( 10 1 rapidly 
becomes rotationally invariant close to the origin. In 
terms of the new variable ^ — (,/2y/a, it becomes 



(11) 



In this limit, we can derive a closed expression for the gap 
probability . Because of the rotational symmetry we 
choose J to be a semi-circle of radius r = |^| and obtain 
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where we have introduced the incomplete Bessel function 

/lfl(a;) = ELo(2;/2)^"+''/'^!(" + v)\ for ^ > 0, and zero 
otherwise. Our quenched expression Eq. (12 1 general- 



izes the corresponding result of Ref. [1] for the non-chiral 
Ginibre ensemble, which is given in terms of incomplete 
exponentials e(,[x) = X^n^o-^"/"'- 

Denoting each factor in Eq. ( 12 1 by 1 — A^, expressions 
for the Ek{r) easily follow in terms of the Xe [30]. The ra- 
dially ordered eigenvalue distributions are then obtained 
from the Ek{r) via Eq. ([s]), leading to 
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Figure [2] shows that the individual eigenvalue distribu- 
tions Pkir) nicely add up to the density Eq. ( 11 ). 
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FIG. 2: Quenched spectral density Eq. (Ill and distributions 
of the first eight eigenvalues Eq. ( |13[ ), as well as their sum, all 
in the large-a limit, for u = (left) and i/ = 1 (right). 



Comparison with lattice data. We now come to 
the comparison of our analytical results to quenched lat- 
tice QCD data. For details of the simulation we refer 
to Ref. [21]. The gauge fields were generated in the 
quenched approximation on a 4"* lattice at P — 5.1 (see 
[21] for an explanation of these choices) . The Dirac oper- 
ator introduced in Ref. [25 is a generalization of the over- 
lap Dirac operator [31] to /i 7^ 0. It satisfies a Ginsparg- 
Wilson relation [32] and has exact zero modes at finite 
lattice spacing. We can therefore test our predictions 
in different sectors of topological charge v. In Ref. [21] 
complete spectra of the generalized overlap operator were 
computed for several values of /x and large numbers of 
configurations, and these data are used in the compar- 
isons to the RMT results below. We also used the fit 
parameters S and F from Ref. [21] to determine a and 
^, i.e., no additional fits were performed. 

For the contours dJ we again choose semi-circles, for 
all values of a. Since we prefer to show 2D plots we 
have integrated over the phase of the complex number 

= i?e'^ and display only the radial dependence. Results 
for 1/ — 0,1,2 are shown in Fig.[3]for fi = 0.1 and fi = 0.2, 
corresponding to a — 0.174 and a = 0.615, and in Fig.|4] 
for fi = 0.3 and 11 = 1.0, corresponding to a = 1.42 and 
a = 4.51. (The lattice spacing a has been set to unity.) 
For all values of /i we compare the data to the expansion 
Eq. ([6|, in which only the first three terms were used. 




FIG. 3: Integrated distribution Pi{R) = dO Rpi{R,e) of 
the first eigenvalue for = 0, 1, 2 as a function of the radius R 
for /i = 0.1 (left) and /i — 0.2 (right). The solid lines are the 
RMT results from Eq. (|6|, the histograms are the quenched 
lattice data of Ref. [52]. The bending- up of the RMT curves 
for large R is an artifact of using only the first three terms in 
the expansion (|6]), see text. 
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FIG. 4; Same as Fig. [s] but for /i = 0.3 (left) and /i = 1.0 
(right). For fi = 1.0 we also show the exact RMT result in 
the large-a hmit from Eq. (|13[). 



For /i = 1.0 the data were found to be approximately 
rotationally invariant, and we also compare them to the 



exact result in the large-a limit from Eqs. ( 12 1 and ( 13 1 



(Because of the rotational invariance, only the ratio S jF 
could be determined for /i = 1.0 in Ref. [24], see Eq. (111. 
In this case the value of a used in Eq. ([6]) is an extrapo- 
lation, assuming that S is independent of /i.) The agree- 
ment between data and analytical curves is excellent ex- 
cept for 1/ = 1,2 at /Lt = 1.0 (see Fig. |4]). In these two 
cases we have left the range of validity of RMT. 

We emphasize that while the rise of the distributions 
from zero was in principle already tested in Ref. |24j 
through the density (see Fig. [l] or [2]), their decrease repre- 
sents a new, parameter-free test. Note also that because 
of the integration over the phase, the signal is much bet- 
ter than in Ref. [H]. This allows us, for the first time, to 
successfully test the RMT predictions for v = 2. 

Figures |3] and |4] also show the effect of truncating the 
Fredholm expansion ([6]): The analytical curves bend up 
for large R after (almost) touching zero. Higher-order 
terms in the expansion ([6| only affect the tail of the dis- 
tributions. They will "repair" the bending-up and ensure 
that the tails remain zero, just as the data. The same ef- 
fect was observed earlier for real eigenvalue distributions 
P7] . This feature of our approximation can be seen most 
clearly when comparing to the exact result in the large-a 
limit, see Fig. |4] (right), in which we can observe how the 
expansion converges in the case of large a. 

Conclusions. We have shown that the distributions 
of individual complex eigenvalues from non-Hermitian 
RMT agree very well with the corresponding distribu- 
tions of the complex eigenvalues of the quenched QCD 
Dirac operator closest to the origin in three different 
topological sectors. As in the Hermitian case, these 
distributions are much easier to compare with than the 
density, in which a plateau may not be observable due 
to appreciable finite-volume corrections. Our analytical 
results are also relevant for other non-Hermitian systems 
in the chiral unitary symmetry class. In the future, it 
would be interesting to compute (and apply) similar re- 
sults for the orthogonal and symplectic symmetry classes. 
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